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1 Introduction 

It has been shown that for a dynamical relaxation process in which a system, 
evolving according to a dynamics of model A 0, and is quenched from very high 
temperature to the critical temperature, there emerges a universal dynamical scaling 
behaviour already within the short-time regime. This rather unexpected scaling 
seems to exist, since spatial correlations and correlation time diverge simultaneously 
as the system approaches long- time regime at the critical temperature. For the 
short-time regime, the finite-size scaling form of the time evolution of a fc-th moment 
of the magnetization is written as |l[ 

M (k \t,e,m ) = b~ k(3/u M^{t/T{L),b 1/u e,b Xa m ). (1) 

Here b is the scale change, L is the linear dimension of the system, (3, v are the well 
known static critical exponents, r is the autocorrelation time, and e = (T — T c )/T c 
is the reduced temperature. Short-time dynamic behaviour also requires a new 
independent critical exponent xo which is the scaling dimension of the initial mag- 
netization m . It has been shown numerically that dynamic scaling exists || even 
at the very early stages of the relaxation process. 

Rigorous formulation of the finite-size scaling for first-order phase transitions [[|, [| 
resulted in better understanding of the dynamics of the first-order phase transi- 
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tions. In this formalism it has been shown that the phase transition is governed by 
the surface tension between the ordered and disordered phases. The system tun- 
nels between these two metastable states and these transitions are observed during 
simulation studies of the long-term behaviour of the system. For finite systems, 
undergoing first-order transitions, the autocorrelation time r for the relaxation pro- 
cess has been calculated |6|, [7|] for cluster algorithms || U and is given as 

r = L^exp^Z/" 1 ), (2) 

where d is the dimensionality of the system. This form of r can be used to identify 
the order of the phase transition. 



In a series of previous work |10|, pT||, behavioral differences between first- and second- 
order phase transitions have been studied. In these works, empirically distinct 
change in the time evolution of the operator in initial stages of the simulation gave 
a clear indication that first- and second-order phase transitions are grouped into 
two different evolutionary processes. Since short-time dynamic behaviour of second- 
order phase transitions are well understood in terms of dynamic scaling formalism [|TJ , 
in analogy with second-order transition, scaling for first-order transitions may be put 
into more rigorous footing. The success of the finite size scaling arguments and ex- 
plicit form of t given in Eq.(2) led us to study the existence of short-time dynamic 
scaling in first-order phase transitions. In first-order phase transitions the singular- 



ities are governed by the volume of the system. Hence in first-order transitions the 
thermal and magnetic critical indices are replaced by the dimension of the system. 
Combining this information with Eq.(2), we have formulated dynamic scaling form 
of various operators in anology with Eq. (1). In this work, our aim is to show that a 
system exhibiting first-order phase transition obeys well defined dynamic finite-size 
scaling rules during quenching from disordered state to the infinite lattice transition 
temperature. We have studied the short-time relaxation processes by using q = 6 
and 7 state Potts models, which are known to exhibit first-order phase transition. 
In this model we studied the time evolutions of the order parameter, the largest 



cluster and the Binder cumulant 12 . 



2 Model and Method 

The Hamiltonian of the g-state Potts model [0, [IJJ is given as 



-m = E ( 3 ) 

<ij> 

where the spin s can take values 1, . . .,q, f3 = 1/ksT is the inverse temperature, 
K = J I (JcbT), 6 is the Kroneker delta function, and sum is over all nearest-neighbour 
pairs on two-dimensional lattice. In equilibrium the g-state Potts model is exactly 
solvable. The critical point locates at K c = log(l + y/q). In principle, any type 



of dynamics can be given to the system to study non-equilibrium time evolution 
process. In this work, we used nonconserved dynamics of Model A 0. In order 
to study dynamic scaling in systems exhibiting first-order phase transitions the 
following operators are considered: 

1. Moments of the order parameter (M) 



M« = (( q -^-^) k ) (4) 

p a = N a /L d , N a being the number of spins with s = a, L the linear size and 
d is the dimensionality of the system. 

2. Binder cumulant (B) 



B = 1-— — ? (5) 

3M( 2 ) 2 



3. Largest cluster (C„ 



C m = ^- d (N Cm ) (6) 

Nc m is the number of spins belonging to the largest cluster in each configuration. 
Largest cluster gives the time evolution of the average of the largest cluster found in 
each configuration. This quantity scales like the susceptibility. Hence in a first-order 
phase transition, it grows like volume. 

4 



For the first-order phase transitions, since the static critical exponents are replaced 
by the dimension of the system, rather than calculating the static critical indices, 
one can test the validity of the dynamic scaling assumption at the initial stages of 
the simulation and obtain the surface tension as the result of the scaling. For the 
computational simplicity, the initial magnetization mo is set to zero. For second- 
order phase transitions, the finite-size behaviour of the magnetization is given by 
Eq.(l). Here, (3/v = Y H — d. Since Y H and Y T are equal to the dimension of 
the system, for first-order phase transitions, the order parameter (M^), Binder 
cumulant (B) and the largest cluster (C m ) scale according to 

f^t/riLxlO,^) = f L2 (t/r(L 2 ),0,L 2 ) (7) 

where r(L) is autocorrelation time of the lattice with the linear size L. Application 
of this form to data can show scaling for various size lattices. In the following section 
we have presented our results. 

3 Results and Discussions 

Following the considerations started in previous section, we have studied the two- 
dimensional q = 6 and 7 state Potts models evolving in time according to dynamics 
of model A ||. Our main objective is to observe the dynamic scaling, hence we have 
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prepared lattices with vanishing order parameter, avoiding the complications due to 
having an extra parameter xq. This is achived for q = 6 and 7 state Potts models 
by choosing the lattice sizes as the integer multiples of q. Totally random initial 
configurations are quenched at the corresponding infinite lattice transition temper- 
ature. Simulations are performed on 6 different lattice sizes by using Wolff cluster 
update algorithm. For each q and L the averages are taken over 10000 different 
samples. Errors are calculated by dividing the samples into ten subsamples. As the 
lattice size grows, number of iterations for thermalization grows according to grow- 
ing tunneling time (Eq. 2). For q = 7 and larger lattices up to 30000 iterations are 
necessary for thermalization. The chosen lattice sizes are L = 42, 60, 72, 90, 96, 102 
and L = 35, 49, 63, 77, 91, 105 for q = 6 and q = 7 respectively. 



The two-dimensional g-state Potts model is known to undergo first-order phase 
transition for q > 4 []nj. Even though the q = 7 state Potts model exhibits 
strong first-order behaviour, the correlation length is about 50 lattice sites. Hence 
for q = 7 [HJ, the largest lattices are expected to show good scaling behaviour 



without any need to correction to scaling terms. For smaller lattices, however, one 
needs to consider the correction to scaling according to the finite-size scaling theory 
for first-order phase transitions. The general form of the corrections to the scaling 
can be given as polynomial in j^, which can be written as 
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(A) L = 4,(1 + ^ + ^ + ...). (8) 

This form indicates that all of the observables scale if one calculates Aq by fitting 
the correction to scaling terms j|, [5], PH . The correction to scaling plays even more 
profound role for q = 6 state Potts model where the correlation length is larger than 
even the largest lattice. The correction to scaling for each observable is obtained by 
fitting the averages, taken over 10000 iterations after the thermalization, to Eq. (8) 
and the expansion coefficients Ai,A 2 ,.. are calculated for q = 6 and 7. 



In Figure l.a, the time evolution of the order parameter is plotted for q = 6. As 
one can observe, for each lattice size, starting from totally random configuration 
m = 0, the order parameter evolves to a plato. For large enough lattices, since 
Yh — d vanishes for first-order phase transitions, one can expect the same long-term 
behaviour for all different lattice sizes. In fact this is the case, within the errorbars, 
for our largest two lattices. In order to see scaling for smaller lattice sizes we have 
performed long runs, after thermalization, and the correction to scaling terms (Eq. 
8) are fitted to the order parameter values. In figures l.a and l.b the raw data and 
the scale form is presented for q = 6 and 7 respectively. Figures l.c and l.d show the 
scaled forms of the data in figures l.a and l.b respectively. For q = 7, the correction 
to scaling is almost negligible for lattices larger than L = 65. 
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M 


0.0084 


± 


0.0011 


0.0185 


± 


0.0015 


B 


0.0076 


± 


0.0004 


0.015 


± 


0.002 


Cm 


0.0084 


± 


0.0011 


0.0175 


± 


0.0015 



Table 1 : 2 a od for q = 6 and 7 



Similarly, for the averages of the maximum cluster, which is expected to grow like 
the volume, we have observed that similar scaling behaviour exists. Figures 2. a and 
2.b are plots of the Monte Carlo data for q = 6 and q = 7 respectively. Figures 2.c 
and 2.d are the scaled form of the above mentioned data. 

The last quantity that we have observed is the Binder cumulant. Binder cumulant 
is a scaling function and also is a ratio of two quantities of equal anomolous di- 
mensionality. Hence, correction to scaling terms are almost negligible even for very 
small lattices. 

These scaling studies enable us to calculate the order-disorder surface tension 2a od . 
The surface tensions 2<j od of q = 6 and q = 7 state Potts models are calculated from 
the autocorrelation time of the relaxation processes for the observables. In Table 
1. we have presented 2<r od which are obtained from the relaxation of three different 
quantities. Depending on the quantity, the value of the surface tension is observed 
to vary slightly. Nevertheless, the surface tension, within errorbars, is 0.008 ± 0.001 
and 0.017 ± 0.004 for q = 6 and q = 7 respectively. The error on the surface tension 
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can be taken as the fluctuation of the values obtained using different operators. 



4 Conclusions 

In conclusion we have numerically simulated the dynamic relaxation process of the 
two-dimensional q = 6 and 7 state Potts models starting from random initial states 
with vanishing initial order parameter. Here in this preliminary work we have inves- 
tigated the dynamical scaling properties of the first-order phase transitions. This 
work is based on two well established facts that the autocorrelation time of the 
critical relaxations in first-order phase transitions are given by the instanton calcu- 
lations J7| and all infinities of the thermodynamic quantities are governed by the 
volume of the system [|], [5], [T7| . Under these assumptions one may expect that any 



thermodynamical quantity exhibits dynamical scaling considering the correction to 
scaling terms. 



We have demonstrated that for first-order phase transitions a universal scaling be- 
haviour emerges already in the macroscopic short-time regime of the dynamical 
evolution. This scaling behaviour resembles closely dynamic scaling which seems to 
exist in second-order phase transitions. Furthermore, such a scaling opens new and 
alternative methods of calculating surface tension and it can be used to distinguish 
weak-first-order phase transitions from the second-order one. 
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Figure captions 

Figure I. (a) and (b) are the time evolution of the order parameter M, and (c) and 
(d) are their scaled form for q = 6 and 7-state Potts model respectively. (The 
errorbars are omitted from the scaled forms for clarity of the figures.) 

Figure 2. Same as fig. 1 but plots are for the maximum clusters. 
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